Sequential deconfinement of quark flavors in neutron stars 
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We suggest a scenario where the three light quark flavors are sequentially deconfined under in- 
creasing pressure in cold asymmetric nuclear matter as found, e.g., in neutron stars. The basis 
for our analysis is a chiral quark matter model of Nambu-Jona-Lasinio (NJL) type with diquark 
pairing in the spin-1 single flavor (CSL), spin-0 two flavor (2SC) and three flavor (CFL) channels. 
We flnd that nucleon dissociation sets in at about the saturation density, no, when the down-quark 
Fermi sea is populated (d-quark dripline) due to the flavor asymmetry induced by /3-equilibrium 
and charge neutrality. At about 3no u-quarks appear and a two-flavor color superconducting (2SC) 
phase is formed. The s-quark Fermi sea is populated only at still higher baryon density, when the 
quark chemical potential is of the order of the dynamically generated strange quark mass. We 
construct two different hybrid equations of state (EoS) using the Dirac-Brueckner Hartree-Fock 
(DBHF) approach and the EoS by Shen et al. in the nuclear matter sector. The corresponding hy- 
brid star sequences have maximum masses of, respectively, 2.1 and 2.0 Mq. Two- and three-flavor 
quark-matter phases exist only in gravitationally unstable hybrid star solutions in the DBHF case, 
while the Shen-based EoS produce stable configurations with a 2SC phase-component in the core 
of massive stars. Nucleon dissociation via d-quark drip could act as a deep crustal heating process, 
which apparently is required to explain superbusts and cooling of X-ray transients. 

PACS numbers: 12.38.Lg, 26.60.Kp, 97.60.Jd 



I. INTRODUCTION 



The phenomenology of compact stars is intimately con- 
nected to the EoS of matter at densities well beyond 
the nuclear saturation density, tiq = 0.16 fm~'^. Com- 
pact stars are therefore natural laboratories for the explo- 
ration of baryonic matter under extreme conditions, com- 
plementary to those created in terrestrial experiments 
with atomic nuclei and heavy-ion collisions. Recent re- 
sults derived from observations of compact stars provide 
serious constraints on the nuclear EoS, see [1| and ref- 
erences therein. A stiff EoS at high density is needed 
to explain the high compact-star masses and radii, for 
which there is growing evidence from recent observa- 
tions. A mass of M '^ 2.0 M0 has been reported for 
some low-mass X-ray binaries (LMXBs), e.g., 4U 1636- 
536 (^ , based on the assumption that the abrupt drop in 
the coherence of the lower kilohertz quasiperiodic oscil- 
lation (QPO) may be related to the innermost stable cir- 
cular orbit, see also Q. From observations of the bright 
isolated neutron star RX J1856. 5-3754 (shorthand: RX 
J1856) in the optical and X-ray frequency range, a con- 
servative lower limit of the apparent neutron star radius 
of Rao = 16.5 km is derived [j|. This corresponds to a 



'Electronic address: blaschke@ift.uni .wroc.pT] 
t Electronic address: fredrik.sandin@g niail.com] 
-^Electronic address: thomas@ift.uni. wroc.pl| 
§ Electronic address: p ens . berdermann@desy. de | 



true (de-redshifted) radius of i? = 14 km for a 1.4 M© 
neutron star, or equivalently, to a star mass of at least 
2.1 Mq when the radius does not exceed 12 km [l|. An- 
other example is EXO 0748-676, an LMXB for which 
the compact-star mass and radius have been constrained 
to M > 2.10 ± 0.28 Mq and R > 13.8 ± 0.18 km Q 
by a simultaneous measurement of the Eddington limit, 
the gravitational rcdshift, and the flux of thermal ra- 
diation. However, the status of the results for the lat- 
ter object is unclear, because the gravitational redshift 
z = 0.35 observed in the X-ray burst spectra [y| has 
not been confirmed, despite numerous attempts. Further 
constraints on the masses and radii of compact stars have 
been reported (see , e.g., 0,13), but they deserve a care- 
ful discussion which is beyond the scope of the present 
paper. While compact-star phenomenology apparently 
points towards a stiff EoS at high density, heavy-ion col- 
lision data for kaon production p and elliptic flow [l^l 
set an upper limit on the stiffness of the EoS [11] . 

A key question regarding the structure of matter at 
high density is whether a phase transition to quark mat- 
ter occurs inside compact stars, and whether it is ac- 
companied by unambiguous observable signatures. It 
has been argued that the observation of a compact star 
with high mass and large radius, as reported for EXO 
0748-676, would be incompatible with a quark core 01, 
because quark deconfinement softens the EoS and low- 
ers the maximum mass and corresponding radius. How- 
ever, Alford et al. have demonstrated [ll| with a few 
counter examples that quark matter cannot be excluded 
by this argument. In particular, for a recently developed 



"Id 


(Id 


f. 


m^ 











He 



(Id Jl. Jl 



Hd Jl. J" 



pure nuclear it 
(NM) 



d-quark drip 
(NM + d-CSL phase) 



2-flavor quark ma 
(2SC phase) 



3-navor quark mauer 
(CFL phase) 



FIG. 1: (Color online) Schematic picture of chemical poten- 
tials (columns) and sequential deconfinement of quarks with 
increasing baryon density (from left to right). The flavor de- 
pendent thresholds for chiral symmetry restoration (decon- 
finement) are approximately given by the dynamically gen- 
erated quark masses ruf, f = u,d,s (dashed lines). With 
increasing quark chemical potential, fj, — {jj,u + iid)l'i-, the 
d-quark chemical potential is the first to reach the threshold 
in isospin asymmetric matter. Nucleon dissociation therefore 
sets in as d-quarks are deconfined. Still higher \i is needed to 
form 2-flavor and 3-fiavor quark matter phases. 



hybrid star EoS [I4I, based on the DBHF approach in 
the nuclear sector and a three-flavor chiral quark model 
[lj|, stable hybrid stars with masses ranging from 1.2 
Mq to 2.1 Mq is obtained; in accordance with modern 
mass-radius constraints, see also [IJ]. In this model, a 
sufficiently low critical density for quark deconfinement 
has been achieved by a strong diquark coupling, while 
a repulsive vector meanfield in the quark matter sector 
resulted in sufficient stiffness to achieve a high maximum 
mass of the compact star sequence. The corresponding 
hybrid EoS for symmetric matter was shown to fulfill the 
constraints derived from elliptic flow in heavy-ion colli- 
sions. In the present work we discuss a new scenario 
that comprises a sequential transition from nuclear mat- 
ter to deconfined quark matter, which could play an im- 
portant role in asymmetric matter, in particular for the 
phenomenology of compact stars. 

Chiral quark models of NJL type with dynamical 
chiral-symmetry breaking have the property that the 
symmetry is restored (and quarks are deconfined) sep- 
arately for each flavor. When solving the gap and 
charge-neutrality equations self-consistently, the chiral- 
symmetry restoration for a given flavor occurs when the 
chemical potential of that flavor reaches a critical value 
that is approximately equal to the dynamically gener- 
ated quark mass, [if = ficj ~ "mf, where / = u,d,s. 
In asymmetric matter the quark chemical potentials are 
different. Consequently, the NJL model behavior sug- 
gests that the critical density of deconfinement is fiavor 
dependent, see Fig. [TJ In this scenario the down quark 
flavor is the first to drip out of nucleons when the density 
increases, followed by the up quark fiavor and eventually 
also by strange quarks. This behavior is absent in sim- 
ple and commonly applied bag-model equations of state, 
because they are essentially flavor blind. 

Under the /3-cquilibrium condition in compact stars the 
chemical potentials of quarks and electrons are related by 



/id = [Is and fJ.d = [J-u + [J-e- The mass difference between 
the strange and the light quark flavors rus ^ mu, 'm-d has 
two consequences: 1) the down and strange quark den- 
sities arc different, so charge neutrality requires a finite 
electron density and, consequently, 2) /i^ > /i„. When in- 
creasing the baryochemical potential the d-quark chem- 
ical potential is therefore the first to reach the critical 
value, fic,d ~ TTT-d, where the chiral symmetry gets (ap- 
proximately) restored in a first-order transition and de- 
confined d-quarks appear. Due to the finite value of fig, 
the u-quark chemical potential is still below [i^.u ~ "i„, 
while the s-quark density is zero due to the high s-quark 
mass. A single-flavor d-quark phase therefore forms in 
co-existence with the positively charged nuclear-matter 
medium. 

Why has this interesting scenario been left unnoticed? 
One reason is that bag models, which are commonly used 
to describe quark matter in compact star interiors can- 
not address sequential deconfinement. Another reason is 
that the single-flavor d-quark phase is negatively charged 
and cannot be neutralized with leptons. It was therefore 
disregarded in dynamical approaches like NJL models 
which in practice are used mainly to describe the de- 
confined and "pure" quark matter phase only. In the 
following we discuss the single-flavor phase for the first 
time, under the natural assumption that the neutralizing 
background is nuclear matter. Since nucleons are bound 
states of quarks, a mixed phase of nucleons and free d- 
quarks could naturally arise when nucleonic bound states 
dissociate (Mott effect). 



II. PHASE TRANSITION TO QUARK 
MATTER: NUCLEON DISSOCIATION 

The task to develop a unified description of the phase 
transition from nuclear matter to quark matter on the 
quark level, as a dissociation of three-quark bound states 
into their constituents in the spirit of a Mott transition 
has not yet been solved. Only some aspects have been 
studied within a nonrelativistic potential model [iSl . Il6l | 
and within the NJL model [l3| ■ Here we consider a chem- 
ical equilibrium reaction on the form n + n ^^ p + 3d, 
which results in a mixed phase of nucleons and down 
quarks once the d-quark chemical potential exceeds the 
critical value. This scenario is analogous to the disso- 
ciation of nuclear clusters in the crust of neutron stars 
(neutron dripline) and the effect may therefore be called 
the d-quark dripline. We approximate the quark and nu- 
cleon components as subphases, which are described by 
separate models. 

For the nuclear matter subphase we use two alter- 
natives: 1) The DBHF approach [H [H, [13, [HI, [H 
with the relativistic Bonn A potential, where the nucleon 
self-energies are based on a T-matrix obtained from the 
Bethe-Salpeter equation in the ladder approximation. 2) 
The EoS by Shen et al. [2^ , which is based on relativistic 
mean-field theory and includes the contribution of heavy 



nuclei, described within the Thomas-Fermi approxima- 
tion. Despite its drawbacks this EoS is instructive since 
it is available for a large enough range of densities, tem- 
peratures and isospin asymmetries that it qualifies for 
applications in studies of supernova collapse and pro- 
toneutron star evolution. Only very recently, significant 
progress could be made, e.g., in generalizing the nuclear 
statistical equilibrium approach [24| and in implementing 
a quantum statistical description for cluster formation 
and dissociation (Mott effect) [23. The quark matter 



J 



phase is described within a three-flavor NJL-type model, 
which includes diquark pairing channels [ij, [20, 123, [2a] . 
This approach is justified since the /i > domain of the 
QCD phase diagram is rather poorly understood. A more 
fundamental approach, like solving the in-medium QCD 
Schwin ger- Dyson equations in a concrete QCD model 
[29I . ISOl . l3ll | is demanding and is therefore beyond scope of 
this work. The path-integral representation of the NJL 
partition function is given by 
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where fi = i^s + diag^(|, -i, -^)fJ.Q + Aa^a 4- As/^s is 
the diagonal quark chemical potential matrix and rh = 
diagf (ttIu, md, rrig) is the current-quark mass matrix. For 

a = 0, tq = (2/3)^/^1/, otherwise Tq and Aq are Gell- 
Mann matrices acting in, respectively, flavor and color 
spaces. C = 17^7*^ is the charge conjugation operator and 
q — 9^7'^. The scalar quark-antiquark current-current in- 
teraction is given explicitely and has coupling strength 
Gs- The 3-momentum cutoff, A, is fixed by low-energy 
QCD phenomenology (see table I of [321 )■ The spin-0 
and spin-1 diquark currents are Jdo,a = q^ iC^^TA^Aq 
and JDi = q^iCl-fiXj + 72A5 + 73A2)g. While the rela- 
tive coupling strengths tjuq and rjjyi are essentially free 
parameters, we restrict the discussion to the Fierz values, 
Vdo ~ 3/4 and rjoi ~ 3/8, see [S^- Color superconduct- 
ing phases in QCD with one flavor were first discussed 
in Refs. [3J, [Sa, [3^, where it was also pointed out that 
the gap is of order 1 MeV in the spin-1 color-spin-locking 
(CSL) phase. This feature of the CSL phase is robust. 
See [37| for an analysis of our isotropic ansatz for the spin- 
1 diquark current, |38j for its generalization to the non- 
local case, and [33 for a self-consistent Dyson-Schwinger 
approach. 

The gaps and the renormalized masses are determined 
by minimization of the mean-field thermodynamic po- 
tential under the constraints of charge-neutrality and 
/3-e quili briuin, see Fig. [2] For further details, see 
[l3l . [20, [23, [23] ■ Different phases are characterized by 
different values of the order parameters (masses, gaps, 
etc.) and correspond to different local minima of the 
thermodynamic potential. For a particular choice of the 
baryon chemical potential there may be several local min- 
ima of the thermodynamic potential. The physical solu- 
tion is that with lowest free energy, or. equivalently, high- 



est pressure. In Fig. [3| we show the pressure of various 
phase constructions. Since we use separate models for 
the confined and deconfined states of quarks the dissoci- 
ation of nucleons does not appear automatically within 
the model. Instead, for a given value of the baryon num- 
ber chemical potential, three different phase construc- 
tions are considered: 1) The homogenous charge-neutral 
and /3-equilibrated nuclear matter phase. 2) The ho- 
mogenous charge-neutral and /3-equilibrated quark mat- 
ter phase. 3) A charge-neutral equilibrium mixture of 
nuclear matter and quark matter, or two different quark- 
matter phases. For the models considered here we find 
that the asymmetry in three-flavor quark matter (CFL) 
is so small that it makes little sense to consider inhomoge- 
nous phase constructions. For mixtures of nuclear matter 
with one-flavor (d-CSL) or two-flavor quark matter in the 
2SC or normal quark matter (NQ) phase, however, the 
asymmetry is significantly lower and the pressure higher 
when compared with the homogeneous phases. 

In Fig. |4|we plot the thermodynamically favored phase 
in the plane of baryon and charge chemical potentials. 
The hybrid EoS corresponds to the dash-dotted lines in 
Fig. [4| , i.e., the borders between positively and neg- 
atively charged phases, and they are constructed such 
that the corresponding mixture of nuclear matter, quark 
matter and leptons is charge neutral. At low densities 
a mixture with one-flavor quark matter is favored. At 
higher density, beyond the up-quark threshold, a mix- 
ture of one-flavor and two-flavor quark matter is favored. 
The strange flavor occurs at still higher densities. Note 
that the 2SC phase cannot persist at high j/igl since the 
large difference in the Fermi levels of u and d quarks pre- 
vents their pairing and the two-flavor quark matter is 
therefore in the normal phase (NQ). 
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FIG. 2: (Color online) Solution of the NJL gap equations for 
isospin-asymmetric charge-neutral matter. The upper (lower) 
panel corresponds to the hybrid EoS based on the DBHF 
(Shen) nuclear EoS. The asymmetry at a given value of the 
baryon chemical potential, hb, is different in the two cases, 
because the charge density of nuclear matter depends on the 
model used. 



Using the hybrid EoS we calculate compact-star se- 
quences by solving the Oppenheimer-VolkofF equations 
for hydrostatic equilibrium. The hybrid-star sequences 
fulfill all modern constraints on the mass-radius relation- 
ship, see Fig. \E\ For the DBHF hybrid EoS all stars 
with DBHF-I-CSL matter in the core are stable equilib- 
rium solutions, while the appearance of u-quarks and the 
associated formation of a 2SC subphase renders the se- 
quence unstable. The situation is somewhat different for 
the Shen hybrid EoS, because in addition to Shen+CSL 
stars there are stable solutions with 2SC+CSL matter 
in the core. In both cases configurations with strange 
quarks in the core are unstable. The hybrid star se- 
quences 'masquerade' as neutron stars [40[, because the 
mechanical properties are similar to those of nuclear mat- 
ter stars and the transition from nuclear matter to the 
mixed phase is associated with a relatively small discon- 
tinuity in the density. Unmasking neutron star interiors 
might therefore require observablcs based on transport 
properties, which could be strongly modified in presence 
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FIG. 3: (Color online) Pressure of matter in beta-equilibrium 
for different nuclear-matter models (DBHF - upper panel, 
Shen - lower panel) and phase transition constructions with 
color superconducting quark matter (NJL model). Here 
DBHF (or Shen) -f NJL refers to the mixed phase of nu- 
clear matter and quark matter. At low chemical potential 
the quark-matter phase component is the negatively charged 
d-CSL phase, which lowers the asymmetry of the system and 
thereby gives a higher pressure of the mixed phase. At higher 
chemical potentials (1370 MeV for DBHF-I-NJL and 1230 for 
Shen-|-NJL) there is a transition in the quark sector to a 2SC 
phase component. 



of color superconductivity. It has been suggested to base 
such tests of the structure of matter at high density on an 
analysis of the cooling behavior [4l|, |42, |4^, |4J] or the sta- 
bility of rapidly spinning stars against r- modes [45l . |46| . 
It has turned out that these phenomena are sensitive to 
the details of color superconductivity in quark matter. 

The down-quark chemical potential exceeds that of up 
quarks in asymmetric nuclear matter and, as we have 
illustrated above, this could lead to sequential deconfine- 
ment of the quark flavors. We have checked that the 
breaking of the \J(X)a symmetry with a six-point 't Hooft 
interaction does not rule out the single-flavor d-CSL so- 
lution, but the phase border is shifted to higher |^q|, i.e., 
more asymmetry is needed to realise the d-CSL phase in 
that case. A potential consequence of this is that the 
asymmetry of charge-neutral nuclear matter is less than 
that of d-CSL matter with broken [/(I) a symmetry, and 
that the nuclear phase therefore is thcrmodynamically fa- 
vored. As the origin of the V(X)a anomaly is unknown, 
see, e.g., the discussion in [l3|, and the critical asymme- 
try depends on the parametrization of the NJL model 
and on the nuclear matter model used, a definite answer 
whether the d-CSL phase is realised is a matter of further 
investigation. Other effects of the inhomogenous phase 
mixture, e..g.; Coulomb interactions and surface tension 
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FIG. 4: (Color online) Phase diagrams in the plane of baryon 
and charge chemical potential. The dash-dotted line de- 
note the border between oppositely charged phases. The nu- 
clear matter EoS is DBHF (upper) and Shen et al. (lower). 
Only one-flavor and two-flavor solutions are displayed, be- 
cause three-flavor matter is charge-neutral for fiq ~ 0. The 
transition to the nearly symmetric three-flavor CFL phase 
occurrs at fiB ~ 1600 MeV, see Fig. [5] At low densities the 
mixture of nuclear matter with one-flavor quark matter (d- 
CSL) is favored. At higher densities, beyond the up-quark 
threshold, a mixture with two-flavor quark matter is favored. 
The two-flavor phases considered here are the normal quark 
matter (NQ) and the superconducting (2SC) phase. 



should also be considered in a future detailed investiga- 
tion. Irrespective of these unsettled issues it is clear that 
the free energy of the d-CSL phase decreases with in- 
creasing asymmetry, in direct contrast to the behaviour 
of traditional phases such as the nuclear matter phase. In 
the following we discuss another interesting feature of the 
d-CSL phase, which could have important consequences 
for the phenomenology of compact stars. 



III. BULK VISCOSITY AND URCA 

EMISSIVITY OF THE SINGLE-FLAVOR CSL 

PHASE 

Rotating compact stars would be unstable against r- 
modes in the absence of viscosity [43, 13 ■ Constraints 
on the composition of compact-star interiors can there- 
fore be obtained from observations of millisecond pulsars 
pa . l46j . In such investigations the bulk viscosity is a key 
quantity and constraints on matter phases in neutron 
star interiors can be based on its value. Here we consider 
some relevant aspects of the bulk viscosity for color su- 
perconducting phases, starting with the 2SC phase and 
following the approach described in Ref. [43|. Note that 
the 2SC phase considered in [50| is a three-flavor phase, 
for which the nonleptonic process u + d <^ u + s is the 
dominant contribution. This process is not relevant for 
the 2SC phase discussed here, where the strange quark 
Fermi sea is not occupied. 

The temperature-dependent bulk viscosity for the 
2SC-I-CSL phase has been calculated self-consistcntly in 
[51| and is based on the flavor-changing weak processes 
of electron capture and f3 decay 

u + e^ ^ d + Ve , d ^ u + e~ + 9e . (3) 

It has been shown that the bulk viscosity is related to the 
direct URCA emissivity, which for normal quark matter 
was first calculated by Iwamoto |52l | and can be expressed 
as 
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Here Gp is the weak coupling constant and Que is the an- 
gle between the up-quark and electron momenta, which 
is obtained from momentum conservation in the matrix 
element, see Fig. [6l The triangle of momentum conser- 
vation holds for the late cooling stage, when the tem- 
perature is below 1 MeV and neutrinos are untrapped. 
Trigonometric relations are used to find an analytical ex- 
pression for momentum conservation. To lowest order in 
Ode the result is 



PF,d - PF.u - PF,e 



-^PF. 



(5) 



For small angles dde — due, so it is possible to obtain 
an expression for the matrix element of the direct URCA 
process. Following Iwamoto [S^l one has to account either 
for quark-quark interactions to lowest order in the strong 
coupling constant, a^, ^ or the effect of finite masses 



(6) 
i = u,d,e . (7) 
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From ([S])-© and the /3-equilibrium condition, /^^ = /i„ + 
/ie, the angle Ode that determines the emissivity ^ is 
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FIG. 5: (Color online) Compact star sequences (a) and hybrid equations of state (b). The phase structure at the center, r = 0, 
changes with increasing density, as indicated in the figures. Constraints on the compact-star mass come from 4U 1636 [3 and 
on the mass-radius relation from RX J1856 [j. 
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FIG. 6: (Color online) Direct Urea process in quark matter (a) and triangle of momentum conservation for it (b). 



obtained ([5]) 



IV. MECHANISM FOR DEEP CRUSTAL 
HEATING 
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If interactions and masses are neglected, or the Fermi sea 
of one species is closed as in the single-flavor CSL phase, 
it follows that the triangle of momentum concervation in 
Fig. [S] degenerates to a line or can not be closed. In that 
case the matrix element vanishes with the consequence 
that the direct URCA process does not occur, and also 
the bulk viscosity is zero. However, in the mixed nu- 
clear -l-CSL phase there could be important friction and 
pair-breaking/formation processes, which we have not 
yet studied in detail. This could be an interesting is- 
sue for further investigation due to the large difference in 
the masses of baryons and dcconfincd quarks. 



Superbursts are rare, puzzling phenomena observed as 
a extremely long (4-14 hours) and energetic {^ lO'^^erg) 
type-I X-ray bursts from LMXBs. They take place if the 
accreted hydrogen and helium at the surface burns in an 
unstable manner, which is the normal case [53|. As sug- 
gested in [SJ] superbursts could originate from accreting 
strange stars with a thin crust and a core of three-flavor 
quark matter in the color-flavor-locked (CFL) phase. The 
suppression of the neutrino emissivity and heat conduc- 
tivity in the CFL phase [5^, [5^, [S^ , caused by pairing 
gaps that affect all flavors, is of particular importance 
in this supcrburst scenario. Following Cumming et al. 
[58j the underlying mechanism is unstable thermonuclear 
burning of carbon in the crust, at column depths of about 
(0.5 — 3) X 10^^ g cm~^. Carbon is a remnant of ac- 
creted hydrogen and helium at the surface. Observed 
superburst light curves suggest that the burning takes 
place at a depth where the crust reaches temperatures of 
about 6 X 10^ K and column depths of about 10^^ g cm~^. 
Such high temperatures in the crust at a certain depth 
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FIG. 7: (Color online) Density profiles of two stars with 
masses 1.4 M0 and 2.0 M©. In the model adopted here the 
mixed phase of d-CSL quark matter with nuclear matter ex- 
tends up to the crust-core boundary. 



are caused by deep crustal heating [59|, |60|, \&m- The 
important ing redients for the strange-star model of su- 
perbursts [54| are a thin baryonic crust of thickness 100 
to 400 m, an energy release of 1 to 100 MeV per accreted 
nucleon by conversion into strange matter, a suppression 
of the fast direct URCA neutrino emissivity to the order 
of 10^^ erg cm~'^ s~^, and a thermal conductivity, k, of 
quark matter in the range 10^^ — 10^^ erg cm"-'^ s~^ K"-'^. 

In Fig. |7| we show that d-CSL quark matter (in the 
mixed phase with nuclear matter) extends up to the 
crust-core boundary, as strange quark matter does in the 
case of strange stars. One of the main arguments for 
strange matter in the context of a superburst mechanism 
is the fact that superconducting phases, like the CFL 
phase, can suppress fast neutrino emission processes of 
all quark flavors and are able to fulfill the fusion ignition 
condition. This is the case also for the single-flavor CSL 
phase. As we have shown above the fast direct URCA 
process is not possible at all in this phase, while slow 
neutrino cooling processes like bremsstrahlung of elec- 
trons and d-quarks exist. 

We want to estimate the order of magnitude of the en- 
ergy release AE due to partial conversion of ordinary nu- 
clear matter to DBHF-I-CSL hybrid matter at the crust- 
core boundary. As we apply the Gibbs construction of a 
phase transition a density dependent volume fraction of 
the d-quark admixture in the nuclear -I- CSL phase re- 
sults, varying from zero to unity. The caveat of this con- 
struction is that all thermodynamical quantities at the 
onset of the phase transition vary continuously. However, 
in reality an infinitesimally small fraction of the d-quark 
subphase would imply that large residual color forces be- 



tween d-quarks should occur. Therefore, a solution of 
the phase admixture problem with a flnite jump of the 
d-quark admixture at the onset of the tranistion should 
be energetically favored. At the present stage of our work 
we cannot quantify this statement due to the absence of 
conflning forces between color charges in our quark mat- 
ter model. An estimate which we would suggest here is to 
determine the fraction x of dCSL matter at the d-quark 
dripline in the vicinity of the crust-core boundary. Then 
one multiplies the change in energy per baryon due to 
the process n — > ddu with x as an estimate for the prob- 
ablility of this process to occur per accreted nucleon. A 
rough estimate (see Fig. [7]) gives 0.001 < x ^ 0-01 which 
for a jump of the d-quark mass gap by 300 MeV (see 
Fig. [2]) at the chiral transition (d-quark dripline) results 
in 0.6 < AE [MeV]< 6 MeV. This meets well the esti- 
mated range A£' - 1 - 100 MeV [13, ^ and could thus, 
in principle, explain burst ignition at appropriate depths 
for a suitable value of n. 

Therefore, a strange matter core is not necessar- 
ily needed to resolve the superburst puzzle, because a 
hybrid-star model with quark matter in the d-CSL phase 
could have similar properties. Stejner et al. |53| show 
that deep crustal heating mechanisms at the crust-core 
boundary, e.g., conversion of baryonic matter to strange 
quark matter, which can fulflll the constraints of the su- 
perburst scenario provide a consistent explanation also 
for the cooling of soft X-ray transients. Along the lines 
of this argument, we suggest that the d-quark drip effect, 
which leads to a mixture of nuclear matter with single- 
flavor quark matter in the CSL phase, can serve as a deep 
crustal heating mechanism. Superbursts and the cooling 
of X-ray transients are not only consistent with quark 
matter in compact stars but may qualify as a signature 
for its occurrence! 



V. CONCLUSIONS 

In this paper we suggest a new quark-nuclear hybrid 
model for compact star applications that fulfills modern 
constraints from observations of compact stars. Due to 
isospin asymmetry, down-quarks may "drip out" from 
nucleons and form a single-flavor color superconducting 
(CSL) phase that is mixed with nuclear matter already 
at the crust-core boundary in compact stars. The CSL 
phase has interesting cooling and transport properties 
that are in accordance with constraints from the thermal 
and rotational evolution of compact stars [5ll | . It remains 
to be investigated whether this new compact star com- 
position could lead to unambiguous observational con- 
sequences, and whether it is thermodynamically favored 
also when effects like Coulomb screening and surface ten- 
sion are accounted for. We conjecture that the d-quark 
drip may serve as an effective deep crustal heating mech- 
anism for the explanation of the puzzling superburst phe- 
nomenon and the cooling of X-ray transients. 
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